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ABSTRACT 

Estimates of the interstellar deuterium abundance span a wide range of values. Until recently, 
it was customary to adopt the abundance of deuterium measured in the Local Bubble as rep- 
resentative of the local one. Now, it is becoming unclear whether the true local deuterium 
abundance is significantly higher or lower than this value, depending on the interpretation 
given to current data. It is important to deal with the issue of the deuterium variation and see 
whether it challenges our current understanding of the Galaxy evolution. To this aim, we study 
the evolution of deuterium in the framework of successful models for the chemical evolution 
of the Milky Way able to reproduce the majority of the observational constraints for the solar 
neighbourhood and for the Galactic disc. We show that, in the framework of our models, the 
lowest D/H values observed locally cannot be explained in terms of simple astration processes 
occurring during the Galaxy evolution. Indeed, the combination of a mild star formation and a 
continuous infall of unprocessed gas required to fit all the available observational data allows 
only a modest variation of the deuterium abundance from its primordial value. Therefore, we 
suggest that depletion of deuterium on to dust grains is the most likely physical mechanism 
proposed so far to explain the observed dispersion in the local data. 



Key words: 

abundances 



cosmology: observations - Galaxy: abundances - Galaxy: evolution - ISM: 



1 INTRODUCTION 

Deuterium is created during primordial nucleosynthesis (e.g. Wein- 
berg 1977; Boesgaard & Steigman 1985; Olive, Steigman & Walker 
2000) and then destroyed in stars through D (p, 7) 3 He reactions oc- 
curring at relatively low temperatures (of a few 10 5 K). Since it is 
hard to synthesize this loosely bound isotope in significant quanti- 
ties in any astrophysical environment (Epstein, Arnett & Schramm 
1976; Epstein 1977; Prodanovic & Fields 2003), its abundance in 
galaxies is expected to monotonically decrease with time owing to 
stellar cycling. 

Deuterium is the light nuclide produced during Big Bang Nu- 
cleosynthesis (BBN) whose primordial abundance is most tightly 
tied to the baryon-to-photon ratio, 77 (or, equivalently, to the ra- 
tio of the baryon and critical densities at the present time, fib). In 
principle, once its primordial abundance, (D/H) p , is derived from 
measurements in unevolved systems, the value of fib is precisely 
known. However, there is still significant dispersion among the 
abundances derived for a handful of high-redshift, low-metallicity 
absorption-line systems towards quasi-stellar objects (QSOs) with 
reasonably firm deuterium detections. A more precise r\ determina- 
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tion comes from the observations of the cosmic microwave back- 
ground (CMB) power spectrum coupled to Lya forest data (Spergel 
et al. 2003; see, however, Pettini 2006, his fig. 7). The suggested 
value, 7/10 = 6.1^Q 2, where 7710 = 10 10 ?7 = 274fib/i 2 , when 
coupled to predictions from standard BBN theory, leads to a very 
narrow range for the primordial deuterium abundance, (D/H) p ~ 
2.4-2.8 x 10~ 5 (e.g. Coc et al. 2004). 

Relating (D/H) p to the current local abundance of D meets 
with the non-trivial issue of estimating what the true local abun- 
dance of deuterium actually is. In fact, a large scatter is found 
among nearly fifty D/H measurements towards as many Galactic 
lines of sight. Possible explanations (including deuterium deple- 
tion on to dust grains, or recent infall of unprocessed gas; see e.g. 
Pettini 2006, and references therein; Prochaska et al. 2005) lead to 
different interpretations of the data. 

Up to some years ago, models for the chemical evolution of 
the Galactic disc were used to bound the primordial D mass frac- 
tion. Studies of this kind were motivated by the widely differing (by 
an order of magnitude) observational estimates of (D/H) p published 
in the literature, coming from a few high-redshift QSO absorption 
spectra. It was crucial then to infer (D/H) p from the D abundance 
observed in the solar system and the interstellar medium (ISM) as- 
suming the astration factor required to reproduce the largest set of 
available observational data. In general, it was found that (D/H) p 
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Table 1. D/H measurements towards QSOs. A weighted mean of these data (except measurements towards Q 0347— 3819 and PKS 1937—1009 at z abs = 
3.256; see text for details) gives (D/H) = 2.6 ± 0.3 (the error in the mean is the dispersion about the mean divided by the square root of the number of data 
points), a value consistent with the primordial abundance inferred from the CMB power spectrum and the standard BBN theory. 



QSO Target" logJV(HI) z abs D/H ±lcr [Si/H] [O/H] Reference 



PKS 1937-1009 (I) 


LL 


17.86 ± 0.02 


3.572 


(3.3 ±0.3) Xl0~ 5 


-2.7, -l.9 b 




Buries & Tytler 1998a 


Q 1009+2956 


LL 


17.39 ±0.06 


2.504 


3.98±{5;t? xiCr 5 


-2.4, -2.7 C 




Buries & Tytler 1998b 


HS 0105+1619 


LL 


19.422 ± 0.009 


2.536 


(2.54 ±0.23) Xl0~ 5 


-1.85 


-2.0 


O'Meara et al. 2001 


Q 2206- 199 


DLA 


20.436 ± 0.008 


2.0762 


(1.65 ±0.35) xlO- 5 


-2.23 




Pettini & Bowen 2001 


Q0347-3819 


DLA 


20.626 ± 0.005 


3.025 


(3.75 ±0.25) xlO- 5 


-0.95 ± 0.02 




Levshakov et al. 2002 


Q 1243+3047 


LL 


19.76 ± 0.04 


2.526 


2.42t'o 25 xl0~ 5 




-2.79 ± 0.05 


Kirkman et al. 2003 


PKS 1937-1009 (II) 


LL 


18.25 ±0.02 


3.256 


1.61^x10-5 


-2.0 ± 0.5 




Crighton et al. 2004 


PG 1259+593 


HVC 


19.95 ± 0.05 


0.0 


(2.2 ±0.7) Xl0~ 5 




-0 79+ 12 


Sembach et al. 2004 



a LL = Lyman Limit absorption system; DLA = Damped Lya absorption system; HVC = high-velocity cloud. 
b Blue component, red component. 
c Component 2, component 3. 



values higher than ~4-5 x 10 5 could not be reconciled with 
Galactic chemical evolution (GCE) requirements (Steigman & Tosi 
1992, 1995; Tosi 1996, and references therein; Tosi et al. 1998; 
Prantzos & Ishimaru 2001; Chiappini, Renda & Matteucci 2002), 
with sparse claims of the contrary (e.g. Vangioni-Flam, Olive & 
Prantzos 1994; Scully et al. 1997). 

Recently, Kirkman et al. (2001) obtained Space Telescope 
Imaging Spectrograph (STIS) spectra of the Lya and Lyman Limit 
(LL) regions of the z abs = 0.701 absorption system towards QSO 
PG 1718+4807 (the only example of a QSO absorber with a D/H 
ratio ~10 times the value found towards other QSOs). They found 
that at least part of the extra absorption needed on the blue side of 
the main H I, which was previously ascribed to deuterium, is more 
likely due to contaminating hydrogen (see, however, Crighton et 
al. 2003). Even more recently, a few more low-metallicity, high- 
redshift QSO absorption systems have been searched for deu- 
terium, and all show low D abundances (see next section). This 
piece of evidence, together with the concordance between the rj 
value obtained from the temperature anisotropies in the CMB (as 
well as other astronomical measurements of the power spectrum) 
and that suggested by the standard BBN theory for low (D/H) p val- 
ues, leads to the conclusion that (D/H) p is indeed low, not in excess 
of 3 X 10~ 5 (see, e.g., discussions in Romano et al. 2003; and re- 
cent reviews by Steigman 2004, 2006; and Pettini 2006). At this 
point, it is worth mentioning that an internal tension between the 
abundances of the light elements created in BBN might exist - the 
primordial abundances of 4 He and 7 Li inferred from the observa- 
tions being lower than those expected according to the BBN theory 
constrained by WMAP (see, e.g., fig. 5 of Pettini 2006). While for 
4 He one can appeal to previously unrecognized systematic uncer- 
tainties in the measurements (Olive & Skillman 2004), dealing with 
7 Li in halo stars poses awkward implications to stellar evolutionary 
theories. Yet, the inclusion of non-standard depletion mechanisms 
in stellar evolutionary models promises to reconcile the abundances 
observed in halo dwarfs with the primordial 7 Li abundance inferred 
from the CMB power spectrum and the standard BBN theory (see 
Charbonnel & Primas 2005 for a discussion). 

GCE models are however facing a new problem: connect- 
ing the reasonably well-known, low primordial abundance of deu- 
terium with a presently uncertain local value. In fact, according to 
the newest data, the local abundance of deuterium might be either 
nearly primordial, or significantly lower, quite lower than previ- 



ously assumed [(D/H) L i SM = (1.50 ± 0.10) X 10~ 5 ; Linsky 1998], 
depending on the data interpretation. 

In this paper we will examine these new local D data in the 
framework of Galactic chemical evolution models able to repro- 
duce the vast majority of the available constraints on the disc prop- 
erties. The paper is organized as follows: measurements of deu- 
terium at both high and low redshifts are reviewed in Section 2. 
The dispersion in the local data and its implications are discussed 
in Section 3, where the data are interpreted in the framework of 
complete GCE models. Finally, a discussion is presented in Sec- 
tion 4. 



2 LOCAL AND HIGH-REDSHIFT DEUTERIUM DATA 

Data on high-redshift, low-metallicity Damped Lyman-a (DLA) 
and LL absorption systems set lower limits to the primordial abun- 
dance of deuterium. There is statistically significant scatter in the 
D/H measurements for about half a dozen quasar absorption sys- 
tems with [Si/H] <; —1.0 about the mean (see Table 1). The system 
at ,z a b s = 3.025 towards Q 0347-3819 and that at z abs = 3.256 to- 
wards PKS 1937—1009 are excluded from our analysis because of 
their complex velocity structures, which make them rather unreli- 
able (Steigman 2004, 2006; Pettini 2006). Since the primordial D/H 
is thought to be isotropic and homogeneous and the metallicity of 
the absorption clouds used to measure D/H is low (hence no signif- 
icant D astration through stars should have occurred), the scatter is 
hard to explain. It has been proposed that some early mechanism 
for D astration may be the cause for the reported scatter in D/H 
values (Fields et al. 2001; Prantzos & Ishimaru 2001). Yet, the cor- 
relation between D/H and [Si/H] that one would expect in this case 
is not seen (Fig. 1). On the other hand, unrecognized systematic ef- 
fects may be responsible for an artificial dispersion in the data. In 
this respect, it is interesting to notice that the total neutral hydrogen 
column density N(H I) for the two systems with the highest D/H 
values at [Si/H] —2.0 was measured using the drop in flux at the 
Lyman limit. In all the other systems, the dominant constraint on 
the N(H I) value is from the damping wings of the Lya or Ly/3 lines 
(Pettini & Bowen 2001; Crighton et al. 2004). Still, the placement 
of the continuum in the systems with the highest 7V(H I) might be 
compromised by the presence of interlopers, and lead to erroneous 
N(Hl) determinations (Steigman 1994, 2006). 

Further measurements of deuterium in metal-poor environ- 
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Figure 1. Available measurements of D/H in QSO absorption systems from 
Table 1. The abundance of deuterium is plotted against metallicity (ei- 
ther [Si/H] - open symbols, or [O/H] - filled ones). There are no hints of 
decreasing D/H with increasing metallicity. Rather, a large dispersion is 
present across the whole metallicity range. The weighted mean of the data 
(except for measurements towards Q0347-3819 and PKS 1937-1009 at 
z abs = 3.256; see text for details) is also shown (solid line) with its errors 
(dashed lines). 



merits include that in the high-velocity, low-metallicity gas cloud 
Complex C falling on to the Milky Way, also listed in Table 1. 
Using Far Ultraviolet Spectroscopic Explorer (FUSE) and Hub- 
ble Space Telescope (HST) observations of the QSO PG 1259+593, 
Sembach et al. (2004) detect Dl Lyman series absorption in 
Complex C and derive (D/H) C om P io X c = (2.2 ± 0.7) x 10~ 5 , 
(0/H)compi cx c = (8.0 ± 2.5) x 10~ 5 , and (D/0) C om P icx c = 
0.28 ± 0.12. A value consistent with the primordial deuterium 
abundance, (D/H) = (2.3 ± 0.4) x 10~ 5 (weighted mean for three 
regions at Galactic longitudes I = 171°, 183°, 195°), is found also 
from the 327 MHz D line, for a Galactocentric distance of Rq = 
10 ± 1 kpc (Rogers et al. 2005). 

Precise measurements of D/H in the local interstellar medium 
(LISM) were obtained with Copernicus (e.g. Rogerson & York 
1973), HST (e.g. Linsky 1998), and the Interstellar Medium Ab- 
sorption Profile Spectrograph (IMAPS; Jenkins et al. 1999; Son- 
neborn et al. 2000). In recent years, FUSE has added many deter- 
minations of D/H (as well as D/O and D/N) along several lines of 
sight probing the neutral ISM up to ~2.7 kpc away (e.g. Moos et 
al. 2002; Hebrard & Moos 2003; Wood et al. 2004; Hebrard et al. 
2005; Oliveira et al. 2006). As data accumulated, the situation got 
complicated. Within the Local Bubble (LB; a low-density region 
within a distance of ~100 pc in which the Sun is embedded), D/H 
is nearly constant. However, while Wood et al. (2004) state that 
(D/H) lb = (1.56 ± 0.04) x 10~ 5 , a value derived from N(DI) 
and ./V(HI) measurements towards 16 targets, Hebrard & Moos 
(2003), who rely on measurements of D/O and O/H, find (D/H)lb = 
(1.32 ± 0.08) x 10~ 5 , a value significantly lower. At logAMLII) 
i; 20.5 (i.e., distances greater than 500 pc), the few data points 
available up to now display a low deuterium abundance, D/H = 



(0.85 ± 0.09) x 10" , while in the intermediate range [d ~100- 
500 pc; log N(H I) = 19.2-20.7] the D/H ratio is found to vary from 
~ 0.5 x 10~ 5 to - 2.2 x 10~ 5 (see, e.g., Linsky et al. 2005). This 
behaviour is interpreted as due to deuterium depletion on to dust 
grains (Linsky et al. 2005; see also Prochaska et al. 2005; Oliveira 
et al. 2006), following an original idea by Jura (1982). According 
to this picture, the most representative value for the total (gas plus 
dust) D/H ratio within 1 kpc of the Sun would be (D/H)lism > 
(2.19 ± 0.27) x 10~ 5 (Linsky et al. 2005). On the other hand, on 
the basis of D/O and O/H measurements, Hebrard & Moos (2003) 
suggest that the low D/H value at 'large' distances truly reflects the 
present-epoch D/H. In their scenario, (D/H)lism should be lower 
than 1 x 10~ 5 . 

Solar system observations of 3 He permit an indirect determi- 
nation of the deuterium abundance in the Protosolar Cloud (PSC; 
Geiss & Reeves 1972). Such an estimate of the deuterium abun- 
dance 4.5 Gyr ago, (D/H) PSC = (2.1 ± 0.5) x 10~ 5 (Geiss & 
Gloeckler 1998), is remarkably similar to both the primordial and 
the local D/H values, if the highest D/H values measured by FUSE 
represent the actual abundance of deuterium in the vicinity of the 
Sun. It is worth noticing that measures of deuterium in the atmo- 
sphere of Jupiter using the Galileo Probe Mass Spectrometer give 
a similar value, (D/H) PS c = (2.6 ± 0.7) x 10~ 5 (Mahaffy et al. 

1998) . 

Finally, detection of deuterium in a molecular cloud at 
10 parsecs from the Galactic Centre indicates that (D/H)gc = 
(1.7 ± 0.3) x 10~ 6 (Lubowich et al. 2000), five orders of magni- 
tude larger than the predictions of GCE models with no continuous 
source of deuterium in the bulge (Matteucci, Romano & Molaro 

1999) . This discrepancy is suggestive of some replenishment pro- 
cess, likely infall of gas of cosmological composition (Lubowich et 
al. 2000; Audouze et al. 1976). 



3 GALACTIC EVOLUTION OF DEUTERIUM 

Once the primordial abundance of deuterium is settled thanks to (ei- 
ther direct or indirect) observations, the knowledge of the present- 
day Galactic deuterium abundance is needed to establish the degree 
of astration suffered by deuterium in the Milky Way. This in turn 
has profound implications for our understanding of the mechanisms 
of the formation and evolution of the Galaxy. 



3.1 The solar neighbourhood 

If(D/H) p = (2.6±0.3) x 10" 5 ,and(D/H) L isM > (2. 19 ±0.27) x 
10~ 5 according to Linsky et al. (2005), then a low astration factor, 
/ < 1.2 ± 0.3, is required. On the other hand, if (D/H)lism = 
(0.85 ± 0.09) x 10~ 5 , as argued by Hebrard & Moos (2003), then 
/ ~ 3.1 ± 0.7. Here / is the astration factor by number, 



/ = 



(D/H) p 



(1) 



(D/H)lism ' 

In previous work (Romano et al. 2003), we showed that the as- 
tration factor / ~ 1.5 which is required in order to fulfil all the 
observational constraints available for the Milky Way, well repro- 
duces also the PSC and LISM deuterium abundances, provided that 
(D/H)lism = (1.50 ± 0.10) x 10~ 5 t. Such a low astration factor 

t Often quoted in GCE works are the astration factors by mass, F = 
Xt>, p/Xn, t now ' where X v, p and X ~D, t no „ are the abundances of deu- 
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Figure 2. Theory versus observations. The abundance of deuterium measured in QSO absorbers (circles) is plotted against metallicity (either [Si/H] - upper 
panel, open circles; or [O/H] - lower panel, filled circles). Only the systems with firm D detections are considered (see Table 1 for references and Section 2 
for discussion). Also shown are the deuterium abundance in the PSC (Geiss & Gloeckler 1998; upper and lower panels; Sun symbol) and those representative 
of the true D/H in the LISM (lower panel, squares), according to either Linsky et al. (2005; 'high' value) or Hebrard & Moos (2003; 'low' value). The values 
of O/H in the LISM are those measured by Oliveira et al. (2005) for the neutral ISM inside the LB, log(0/H)+12 = 8.54 ± 0.02, and Esteban et al. (2004) for 
Orion (ionized gas plus dust), log(0/H)+12 = 8.73 ± 0.03. Notice that the Oliveira et al. measurement refers to the gas phase only, and thus misses ~25 per 
cent of the total oxygen due to depletion on to dust grains. Theoretical predictions (upper panel: D/H versus [Si/H]; lower panel: D/H versus [O/H]) refer to 
models for the solar vicinity computed with different prescriptions about the stellar lifetimes and IMF: Scalo's (1986) IMF and Maeder & Meynet's (1989) 
stellar lifetimes (solid line); Kroupa et al.'s (1993) IMF and Maeder & Meynet's (1989) stellar lifetimes (dot-dashed line); Scalo's (1986) IMF and Schaller et 
al.'s (1992) stellar lifetimes (dotted line); Kroupa et al.'s (1993) IMF and Schaller et al.'s (1992) stellar lifetimes (dashed line). See Romano et al. (2005a) for 
details about the models. All ratios are normalized to the solar abundances of Asplund et al. (2005). 



is due to the combination of a moderate star formation and a con- 
tinuous infall of external gas during the whole Galactic disc evo- 
lution. Noticeably, this result is almost independent of the specific 
GCE code used (either that developed by Tosi 1988a,b; or that de- 
veloped by Chiappini et al. 1997, 2002). Higher astration factors 
(/ ~ 2-3) were suggested in the past (e.g. Galli et al. 1995; Prant- 
zos 1996; Fields 1996; but see also Chiappini et al. 1997; Tosi et 
al. 1998). The low one reported here is mainly due to the adoption 
of an updated stellar metallicity distribution, which needs more in- 
fall in order to be reproduced. Having that much infall, it comes 
out really hard to get larger astration factors. As stressed above, in 
Romano et al. (2003) we adopted the mean D/H value measured 
by Linsky (1998) in the LB, (D/H) lb = (1.50 ± 0.10) x 10" 5 , as 
representative of the deuterium abundance in the solar vicinity at 
the present time. This assumption seems now to be no longer valid 
in the light of recent, contrasting interpretations of the dispersion 

terium by mass at t = and t = t n0 w = 13.7 Gyr, respectively. F is larger 
than / by a few per cent, because of hydrogen burning into helium and 
heavier species in the course of the Galaxy's evolution. 



in the local data (see discussions in Section 2). Therefore, in this 
section we recompute the evolution of D/H. 

First, we adopt various prescriptions for the stellar lifetimes 
and initial mass function (IMF) in a successful model for the chem- 
ical evolution of the Milky Way. Details about the model can be 
found in Chiappini et al. (2002) and Romano et al. (2005a), where 
results relevant to several chemical species are discussed. Notice 
that here we consider only those (IMF, stellar lifetimes) combina- 
tions which proved to guarantee a good fit to all the observational 
constraints available for the Milky Way (see Romano et al. 2005a). 
The aim is to associate a 'theoretical error' to our estimate of /. 

In Fig. 2 we show model predictions for D/H versus [Si/H] 
(upper panel) and D/H versus [O/H] (lower panel) in the solar 
neighbourhood obtained with: 

(i) the Scalo (1986) IMF and the Maeder & Meynet (1989) stel- 
lar lifetimes (solid line); 

(ii) the Kroupa et al. (1993) IMF and the Maeder & Meynet 
(1989) stellar lifetimes (dot-dashed line); 

(iii) the Scalo (1986) IMF and the Schaller et al. (1992) stellar 
lifetimes (dotted line); 



© 2006 RAS, MNRAS 000, 1-10 



Deuterium astration 5 



Table 2. Theoretical astration factors for deuterium. Listed are the astration 
factors by number, /, suitable for comparison with the observations. Dif- 
ferent astration factors are predicted when assuming differing formulations 
for the stellar lifetimes and IMF in the GCE code. 



IMF 


Stellar lifetimes 


Astration factors 


Scalo (1986) 


Maeder & Meynet (1989) 


1.52 


Kroupa et al. (1993) 


Maeder & Meynet (1989) 


1.83 


Scalo (1986) 


Schalleret al. (1992) 


1.39 


Kroupa et al. (1993) 


Schaller et al. (1992) 


1.61 



(iv) the Kroupa et al. (1993) IMF and the Schaller et al. (1992) 
stellar lifetimes (dashed line). 

Notice that the dotted lines lay on the dashed ones, but they end 
up with a larger D/H at present, D/H = 1.9 x 10~ 5 rather than 
D/H = 1.6 x 10~ 5 . Indeed, different D/H and O/H ratios are 
reached at late times according to the different IMFs and stellar 
lifetimes adopted. In particular, the higher the mass of massive 
stars (responsible for a quick recycling of the gas) and the lower 
the mass of very-low mass objects (which just lock-up material, 
from the point of view of GCE), the lower (higher) the D/H (O/H) 
ratio predicted at present. The fraction of stars with initial mass 
m > 8 Mq (m < 1 Mq) is larger (smaller) for the Kroupa et 
al. (1993) IMF than for the Scalo (1986) one (see Romano et al. 
2005b, their figure 1). This explains the different behaviours pre- 
dicted by the models shown in Fig. 2. A primordial abundance 
of (D/H) p = 2.6 x 10 -5 is assumed for all the models, i.e. the 
weighted mean of the most reliable QSO absorber data available 
in the literature (see Table 1). This value agrees well with that in- 
ferred from CMB anisotropy measurements and the standard BBN 
theory. Also shown in Fig. 2 are D/H measurements from QSO 
absorber spectra (circles; see Section 2 for references), the PSC 
deuterium abundance after Geiss & Gloeckler (1998; Sun symbol), 
and the representative (D/H)lism value according to either Linsky 
et al. (2005; 'high' value) or Hebrard & Moos (2003; 'low' value). 
Two [O/H] values are displayed for the LISM, namely that mea- 
sured by Oliveira et al. (2005) for the neutral medium inside the 
LB, log(0/H)+12 = 8.54 ± 0.02, and that given by Esteban et al. 
(2004) for Orion, log(0/H)+12 = 8.73 ± 0.03, which represents 
the oxygen content of the ionized matter and dust. Notice that part 
of the local oxygen is depleted on to dust grains (~25 per cent; 
see Oliveira et al. 2005; Meyer, Jura & Cardelli 1998; Andre et al. 
2003; Cartledge et al. 2004), which brings the two measurements 
into agreement. All the ratios are normalized to the solar oxygen 
and silicon abundances from Asplund et al. (2005). An apparently 
striking feature of all the models is that only minor deuterium de- 
struction is predicted as [O/H] varies on more than two orders of 
magnitude, increasing from [O/H] < —3 dex up to [O/H] ~ —0.5 
dex. The deuterium abundance is sensibly reduced afterwards. This 
behaviour is easily explained by the large oxygen yield from the 
first Type II supernovae, which leads to a prompt oxygen enrich- 
ment of the ISM only a few million years after the onset of the 
star formation. Indeed, we find [O/H] > —3.0 dex already 4 mil- 
lion years after the onset of the star formation in the Galaxy, which 
raises to [O/H] = —0.5 dex in the following ~300 million years. In 
the meantime, continuing infall of gas of primordial chemical com- 
position helps to keep the D abundance near its primordial value - 
infall of gas is particularly strong during the early Galaxy evolu- 
tion, so that the consumed D is being replenished fast by infall at 
those earlier times. Later on, when the contribution from low- and 
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Figure 3. Evolution of deuterium with time in the solar neighbourhood, 
as predicted by two distinct codes for the chemical evolution of the Milky 
Way. The codes assume different star formation and infall laws, but the 
same IMF (Kroupa et al. 1993) and stellar lifetimes (Schaller et al. 1992 - 
solid and dashed lines; see text for details). Also shown are the predictions 
from a code assuming Scalo's (1986) IMF and Schaller et al.'s (1992) stel- 
lar lifetimes (dotted line). Data: vertical bar at t = 0: mean (D/H) p from 
QSO absorber observations (see Section 2 for references); Sun symbol: 
PSC D abundance from Geiss & Gloeckler (1998); squares: fiducial 'true' 
(D/H)lism value after Linsky et al. (2005; 'high' value) and Hebrard & 
Moos (2003; 'low' value). 

intermediate-mass stars (which are not net oxygen producers) to the 
chemical evolution of the Galaxy becomes increasingly important 
and the infall rate is reduced, a steeper behaviour on the D/H versus 
O/H plot is observed. 

In Table 2 the astration factors by number for all the models 
displayed in Fig. 2 are listed and compared to each other. All the 
models well reproduce the solar abundances of oxygen and silicon 
at t = t a = 9.2 Gyr. A model in which the Scalo (1986) IMF 
and the Schaller et al. (1992) stellar lifetimes are adopted is able to 
destroy D at a level consistent with that suggested by Linsky et al. 
(2005). None of our models is able to predict astration factors larger 
than / — 2. This leads us to favour the dust depletion scenario. 

Deuterium astration in galaxies is expected to be strongly de- 
pendent on the assumptions for the two competing processes of star 
formation and gas infall. In Fig. 3 we compare the results obtained 
for D/H versus time by two GCE codes, the one of Chiappini et 
al. (2002; dashed line) and that by Tosi (1988a,b; solid line), once 
adopting the same prescriptions on the stellar lifetimes (i.e. Schaller 
et al. 1992) and the stellar IMF (Kroupa et al. 1993). The two mod- 
els assume different prescriptions for the star formation and infall 
rates, yet the predicted deuterium evolution is almost the same. This 
is because what really matters is the interplay between subtraction 
and replenishment of gas at each time, and this interconnection is 
severely constrained by a number of independent observables (see, 
e.g., Matteucci 2004). According to Fig. 3, only a mild deuterium 
destruction is predicted in the last 4.5 Gyr. This is consistent with 
the mild Galactic evolution at late times suggested by the small in- 
crease of the overall metallicity of the gas from the time of Sun for- 
mation up to now (Esteban et al. 2004, and references therein). Also 
the observed G- and K-dwarf metallicity distributions clearly show 
that the majority of the stars in the solar neighbourhood formed 
at low metallicities, [Fe/H] < 0, thus pointing again to a low star 
formation activity (and, hence, less evolution) at late times. Also 
shown in Fig. 3 is the Chiappini et al. (2002) model with Scalo's 
(1986) IMF and Schaller et al.'s (1992) stellar lifetimes, which best 
reproduces both the solar and local 'high' D/H values (dotted line). 
It is worth emphasizing that with these prescriptions the model also 
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Figure 4. Models aimed at reproducing the lowest (D/H)lism values 
disagree with other important observational constraints. For instance, the 
model displayed in this figure (solid lines) assumes a very short time scale 
for thin-disc formation in the solar neighbourhood (rrj = 1.5 Gyr), at vari- 
ance with common assumptions (see Chiappini et al. 1997). Hence, while 
a lower (D/H)lism can be attained (upper panel), the G-dwarf metallic- 
ity distribution of solar neighbourhood stars can not be reproduced any 
more (lower panel). Also shown for comparison are the predictions of a 
successful model (td = 7 Gyr; dotted lines), which well reproduces the 
shape of the observed G-dwarf metallicity distribution. This model leads to 
a higher (D/H)lism value. The observed G-dwarf metallicity distributions 
(histograms in the lower panel) are from Wyse & Gilmore (1995; dashed 
histogram), Rocha-Pinto & Maciel (1996; dotted histogram) and J0rgensen 
(2000; solid histogram). Also shown is the K-dwarf metallicity distribu- 
tion by Kotoneva et al. (2002; big dots). The theoretical distributions are 
convolved with a Gaussian (cr = 0.1) to account for the observational and 
intrinsic scatter. 



well reproduces all the other observational constraints (Romano et 
al. 2005a). 

In general, our models predict (D/H)lism values in the range 
1.4-2.0 x 10~ 5 . We need to enhance the ratio of the star formation 
to gas infall efficiencies if a larger fraction of the initial deuterium 
has to be destroyed. But can a more efficient star formation and/or 
less effective infall of unprocessed gas be accommodated within 
successful GCE models? It is immediately seen that the almost con- 
stant metal abundance of the ISM from the time of Sun's formation 



up to now (Esteban et al. 2004) sharply contrasts with rapid star 
formation and/or absence of external sources of unenriched gas. In 
Fig. 4, we show the outputs of a model especially designed to pro- 
duce a significant drop in deuterium abundance from its PSC value 
to the local one (for further modeling and discussions see Tosi et 
al. 1998). This model, assuming td = 1.5 Gyr rather than 7 Gyr 
for the thin-disc formation time scale, fits the solar and local 'low' 
D/H values (Fig. 4, upper panel, solid line), as well as other obser- 
vational constraints (such as for instance the current gaseous and 
stellar mass surface densities at the solar position). However, it is 
barely consistent with the metal abundances of the Sun, it pred- 
its a sharp increase of the metal abundance of the ISM in the last 
4.5 Gyr (at variance with the observations), and it can not account 
for the metallicity distribution of solar neighbourhood G- and K- 
dwarf stars (Fig. 4, lower panel, solid line). The metallicity distri- 
bution of long-lived stars in the solar neighbourhood is a fundamen- 
tal constraint for GCE models, which mainly constrains the infall 
time scale and strength (see, e.g., Chiappini et al. 1997; Matteucci 
2004). In Fig. 4, lower panel, we show as a dashed, dotted, and solid 
histogram the G-dwarf metallicity distributions observed by Wyse 
& Gilmore (1995), Rocha-Pinto & Maciel (1996) and J0rgensen 
(2000), respectively. Also shown is the K-dwarf metallicity distri- 
bution by Kotoneva et al. (2002; big dots). K dwarfs have lifetimes 
older than the present age of the Galactic disc, and are thus ideal 
stars for investigating the chemical evolution of the disc. G dwarf 
stars instead are sufficiently massive that some of them have begun 
to evolve away from the main sequence, so that corrections must 
be taken into account when determining their space densities and 
metallicity (Kotoneva et al. 2002). These corrections are at least 
partly responsible for the non-negligible differences among differ- 
ent observational G-dwarf metallicity distributions (in particular, 
the low-metallicity tail of the distribution and the exact position of 
the peak may vary significantly according to different authors). In 
conclusion, in the framework of our models we can not explain at 
the same time both the lowest observed (D/H)lism values and the 
observed G- and K-dwarf metallicity distributions. 

3.2 Present-day abundance gradients 

In Fig. 5 we display the Galactic gradients of D/H, D/O, and O/H at 
the present time predicted by adopting the prescriptions of Model 
A of Chiappini et al. (2001) as far as the formation and evolution 
of the Galactic disc are concerned, and by varying the prescriptions 
for the IMF and the stellar lifetimes (solid line: Scalo's 1986 IMF 
and Maeder & Meynet's 1989 stellar lifetimes; dashed line: Kroupa 
et al.'s 1993 IMF and Maeder & Meynet's 1989 stellar lifetimes; 
dotted line: Scalo's 1986 IMF and Schaller et al.'s 1992 stellar life- 
times - see Romano et al. 2005a). A value of (D/H) p = 2.6 x 10~ 5 
is assumed for all the models. Also shown are the relevant obser- 
vations. The filled squares in the upper panel of Fig. 5 stand for the 
representative D/H value in the LISM according to either Linsky 
et al. (2005; 'high' value) or Hebrard & Moos (2003; 'low' value). 
The filled triangle in the same panel refers to the measurement by 
Rogers et al. (2005) for the outer Galactic disk. The filled squares in 
the middle panel of Fig. 5 stand for the representative D/O ratio ac- 
cording to either Hebrard & Moos (2003; lower value) or Linsky et 
al. (2005; higher value). To obtain (D/0)lism in this latter case, we 
divided (D/H)lism from Linsky et al. (2005) by the oxygen abun- 
dance of Orion (Esteban et al. 2004). At first glance, it seems that 
GCE models explain the lowest D/O value better than the high- 
est one, while the opposite is true for D/H (upper panel). How- 
ever, one should recall that the current oxygen abundance displays 
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Figure 5. Present-day distributions of D/H (upper panel), D/O (middle panel) and log(0/H)+12 (lower panel) along the disc of the Galaxy. Predictions from 
models assuming either the Scalo (1986) IMF and the Maeder & Meynet (1989) stellar lifetimes (solid line), or the Kroupa et al. (1993) IMF and the Maeder 
& Meynet (1989) stellar lifetimes (dashed line), or the Scalo (1986) IMF and the Schaller et al. (1992) stellar lifetimes (dotted line) are shown together with 
the relevant observations. Upper panel, filled squares: representative D/H value in the LISM according to Linsky et al. (2005; 'high' value) and Hebrard & 
Moos (2003; 'low' value); filled triangle: deuterium measurement from the 327 MHz D line in the outer disc (Rogers et al. 2005). Middle panel, filled squares: 
representative D/O value in the LISM according to Hebrard & Moos (2003; 'low' value) and Linsky et al. (2005; 'high' value - in this case the oxygen 
abundance measured by Esteban et al. 2004 for Orion has been taken as representative of the total oxygen content of the LISM). Also shown (big star) is the 
D/O ratio towards LSE 44 measured by Friedman et al. (2006; see discussion in the text). Lower panel, small filled circles: compilation of log(0/H)+12 versus 
Rq data from Chiappini et al. (2001); filled square: oxygen abundance in Orion from Esteban et al. (2004). 



a non-negligible spread. Indeed, several indicators - H II regions, 
B-type stars, the neutral medium - show that the local oxygen abun- 
dance varies in the range 8.3 < log(0/H)+12 < 9.0 (Fig. 5, lower 
panel, and Oliveira et al. 2006, their table 12). So, in Fig. 5, mid- 
dle panel, we show also the D/O ratio recently measured by Fried- 
man et al. (2006) for the line of sight towards LSE 44. Along this 
line of sight, (D/H)lse 44 = 2.24±°;g° x 10" 5 , consistent with the 
lower limit suggested by Linsky et al. (2005), and (0/H)lse 44 = 
11-3^3'g x 10~ 4 (that is the highest O/H value measured locally 
by FUSE). Now, within the errors, the model adopting the Scalo 
(1986) IMF and the Schaller et al. (1992) stellar lifetimes is in 
satisfactory agreement with both the observed 'high' (D/H)lism 
and 'low' (D/0)lism, provided that (D/0)lism is 'low' as a conse- 
quence of both a high deuterium and a high oxygen abundance in 
the LISM. 

It is seen that the shape of the gradients slightly changes when 
changing the prescriptions about the IMF and the stellar lifetimes. 
In the Galactocentric distance range Rg = 4-12 kpc the theoretical 
gradients of log(D/H), log(D/0), and log(0/H)+12 are 0.03 ± 0.01 
dex kpc" 1 , 0.08 dex kpc" 1 , and -0.055 ± 0.005 dex kpc" 1 , re- 



spectively. In the outermost regions of the disc the gradient of D/H 
flattens, while those of D/O and O/H steepen, owing to the lower 
and lower amount of gas processed by stars outward in the disc. 

Before concluding that, according to GCE model predictions, 
one shall observe D/H ratios close to the primordial value at the 
largest radii in spiral galaxies, it is worth reminding that the model 
predictions depend on a number of uncertain free parameters. 
Model A of Chiappini et al. (2001) produces a fairly steep O/H 
gradient at the largest radii (Fig. 5, bottom panel), while the ma- 
jority of the data suggest a flatter behaviour. A flatter O/H gradient 
in the outer disc is easily obtained by suppressing the gas density 
threshold which regulates the star formation process during the pre- 
ceding halo phase. The solid lines in Fig. 6 stand for Chiappini et 
al.'s (2001) Model C results. This model is analogous to Model A 
(Fig. 5, solid lines), except that it allows star formation in the halo 
to go on even when the gas density drops below the threshold value 
(~ 4 M Q pc" 2 ). The expected D/H ratio is now ~ 2.5 x 10" 5 at the 
outermost radii, i.e. below the adopted primordial value [(D/H) p = 
2.9 x 10" 5 for this model]. Also shown in Fig. 6 are the predictions 
from the model of Tosi et al. (in preparation), which is an updated 
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Figure 6. Same as Fig. 5, but for Model C of Chiappini et al. (2001 ; solid lines) and for the model of Tosi et al. (in preparation), which is an updated version of 
the model of Tosi (1988a,b) extended to the whole disc (dashed and dotted lines). Here (D/H) p = 2.9 X 10~ 5 (solid and dashed lines) or 2.3 X 10 -5 (dotted 
lines), and not 2.6 X 10~ 5 . 



version of the model by Tosi (1988a,b) adopting the Kroupa et al. 
(1993) IMF and the Schaller et al. (1992) stellar lifetimes, and ex- 
tending to the whole disc [dashed lines: (D/H) p = 2.9 X 1CP 5 ; dot- 
ted lines: (D/H) p = 2.3 x 10~ J ]. This model tends to predict higher 
oxygen abundances, hence lower D/O ratios, across the disc. 

We conclude that, in principle, accurate measurements of both 
D/H and O/H in the outermost regions of spiral galaxies can al- 
low us to discriminate among existing scenarios of galaxy for- 
mation and evolution. In fact, only minor uncertainties related to 
the adopted stellar yields affect the interpretation of the data in 
this case, owing to the well-known nucleosynthetic origin of both 
species. Hence, observational efforts in this direction are highly de- 
sirable. 

In summary, notice that all the GCE models discussed above 
favour minor deuterium variations [Alog(D/H)/Ai?G — 0.02 dex 
kpc" 1 ] in the Galactocentric distance range sampled by current 
observations (Rg — 8-11 kpc), hence favouring the highest 
(D/H)lims determinations as representative of the local deuterium 
abundance. 



4 DISCUSSION AND CONCLUSIONS 

In general, observations of deuterium abundances provide lower 
bounds on its primordial abundance. However, at sufficiently high 



redshifts and low metallicities, one expects the primordial abun- 
dance of deuterium to reveal itself as a plateau. Instead, statistically 
significant scatter of the D/H measurements has been found in sev- 
eral high-redshift QSO absorbers. Before claiming that a disper- 
sion around the mean value is actually present, one has to carefully 
check possible systematics affecting the abundance derivation. It 
has been suggested that in the absorbers with the lowest H I column 
densities, interlopers might indeed contribute to the inferred Dl 
column densities, while in those with the highest H I column den- 
sities, interlopers might affect the wings of the HI lines (Steigman 
2004, 2006, and references therein). Yet, it is worth noticing that the 
placement of the continuum and the fitting of the damping wings 
are closely related for DLAs: if one gets the continuum wrong, it 
is unlikely that he could fit the line profile satisfactorily (M. Pet- 
tini, private communication). In this paper, we assume (D/H) p = 
(2.6 ± 0.3) x 10~ 5 , i.e. the weighted mean of D/H measurements 
towards QSOs (indirect determinations give a consistent value; see 
discussions in Sections 1 and 2). 

Of more concern is the dispersion of local D/H data. In the 
LISM, D/H is found to vary from ~ 0.5 x 10~ 5 to ~ 2.2 x 10~ 5 
(Hebrad et al. 2005; Linsky et al. 2005; Williger et al. 2005; Fried- 
man et al. 2006; Oliveira et al. 2006; to quote only the most recent 
works). Basically, two mechanisms are suggested in the current lit- 
erature which might explain the observed variation (Pettini 2006, 
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and references therein; see also Lemoine et al. 1999; Prochaska et 
al. 2005): 

(i) differing ISM conditions along different sight lines, which 
determine different degrees of deuterium depletion on to dust 
grains; 

(ii) localized infall of unprocessed gas, which modifies the deu- 
terium abundance in the surveyed region while leaving unchanged 
those of more abundant species, such as oxygen. 

In the first case, one should recover the true local D/H in interstellar 
clouds which have been accelerated by supernova shocks, since the 
depletion of refractory elements is much reduced relative to that of 
the normal quiescent ISM there; in this case, the true local abun- 
dance of deuterium would be the highest observed one (Linsky et 
al. 2005, and references therein). In the second case, the true local 
abundance of deuterium would be the lowest observed one. In any 
case, the mean D/H value measured within the Local Bubble should 
not be considered any more as representative of the actual degree 
of astration suffered by deuterium in the solar neighbourhood. By 
assuming that the primordial abundance of deuterium is reasonably 
well known, one can usefully bind the deuterium astration factor 
of the solar vicinity. An astration factor (by number) either as low 
as / < 1.2 ± 0.3 or a factor of 2-3 higher is permitted by the 
observations. 

In this paper, we have computed a number of GCE models 
for the solar vicinity as well as for the whole Galactic disc and 
conclude that: 

(i) Only low astration factors (not in excess of / ~ 1.8) 
are compatible with GCE requirements: our models predict 
(D/H)lism = 1.4-2.0 x 10~ 5 starting from a primordial value of 
(D/H) p = 2.6 x 10~ 5 . This is not surprising: in the absence of 
supernova-driven winds which efficiently remove the metal-rich 
ejecta of dying stars from the Galaxy and/or in the absence of a 
peculiar IMF, the relatively high present-day gas content and low 
metallicity are indicative of modest astration (and, hence, modest 
D destruction; see also Tosi et al. 1998). We emphasize here once 
again that the low astration factors are due to the combination of 
moderate star formation and continuous infall of gas which are 
needed in order to reproduce the available Milky Way data. 

(ii) Small variations in the predicted D astration factor are pro- 
duced by changing the prescriptions on the IMF and the stellar 
lifetimes. In particular, by adopting the Scalo (1986) IMF and the 
Schaller et al. (1992) stellar lifetimes we predict a D astration fac- 
tor which clearly favours the low one suggested by Linsky et al. 
(2005). With this choice for the IMF and the stellar lifetimes, the 
model also gives a good fit to many other observational constraints 
for the Galaxy (Romano et al. 2005a). 

(iii) When the model is forced to reproduce the lowest D/H 
values observed locally, the agreement between model predictions 
and relevant observations is lost, for one or more observables (the 
most relevant observational constraints being - in this context - the 
G-dwarf metallicity distribution, the mild increase of the overall 
metallicity of the ISM from the time of Sun's formation up to now, 
the present-day mass and gas surface densities, the present-day in- 
fall rate; see also Matteucci 2004). 

(iv) In order to attain the lowest D/H values observed in the 
LISM, the smooth, gentle decline which nicely accounts for the 
PSC D abundance, must turn into a steeper behaviour during the 
latest phases of Galaxy's evolution. However, a large gas cycling 
through stars in the last ~5 Gyr is unlikely, in the light of the small 
increase of the global metal abundance from the time of Sun's for- 



mation up to now suggested by several independent indicators (Es- 
teban et al. 2004; their table 15 and references therein). 

(v) In principle, joint observations of deuterium and oxygen 
abundances in the outermost regions of galactic discs can shed light 
on the mechanisms of spiral galaxy formation and evolution. 

In conclusion, we favour a scenario in which D/H in the solar neigh- 
bourhood declines mildly during Galaxy's evolution, due to the 
combined effect of moderate star formation and continuous infall 
of external gas. In this framework, the dispersion in the current lo- 
cal D/H data might be at least partly explained by different degrees 
of dust depletion along different lines of sight. 
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